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Abstract 

Helicity transfer in a shell model of turbulence is investigated. In particular, we study 
the scaling behavior of helicity transfer in a dynamical model of turbulence lacking inversion 
symmetry. We present some phenomenological and numerical support to the idea that Helicity 
becomes -at scale small enough- a passively-advected quantity. 
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"non sono nata per condividere I'odio, ma I'amore" 

1 Introduction 

One of the most intriguing problems in three dimensional fully developed turbulence (FDT) is related 
to the appearance of anomalous scaling laws at high Reynolds numbers, i.e. in the limit when Navier- 
Stokes dynamics is dominated by the non-linear interactions. 

The celebrated 1941 Kolmogorov theory (K41) was able to capture the main phenomenological 
ideas by performing dimensional analysis based on the energy transfer mechanism. Kolmogorov 
postulated that the energy cascade should follow a self-similar and homogeneous process entirely 
dependent on the energy transfer rate, e. This idea, plus the assumption of local isotropy and 
universality at small scales, led to a precise prediction on the statistical properties of the increments 
of turbulent velocity fields: 5v(r) ~ \v(x+r) —v(x)\ ~ (r-e(r)) 1 / 3 . From this the scaling of moments 
of Sv(r), the structure functions, can be determined in terms of the statistics of e(r), i.e. 

S p (r) = ((5v(r)Y) = C p ((e(r)f 3 y/\ (1) 

where C p are constants and the scale r is supposed to be in the inertial range, i.e. much smaller 
than the integral scale and much larger than the viscous dissipation cutoff. If S p (r) ~ r^ p ) and 
(e p (r)) ~ r T ^ then 

C(p)=p/3 + r(p/3). (2) 

In the K41 the e(r) statistic is assumed to be r-independent, that is r(p) = 0, implying ((p) = 
|, Vp. On the other hand, there are many experimental and numerical [J results telling us that 
K41 scenario for homogeneous and isotropic turbulence is quantitatively wrong. Strong intermittent 
bursts in the energy transfer have been observed and non trivial r(p) set of exponents measured. 



Many different authors have focused their attention on the possible role played by hehcity, the 
second global invariant of 3D Navier-Stokes eqs. ||, ||, [|, |9], for determining leading or sub-leading 
scaling properties of correlation functions in the inertial range. 

Recently P, [1(J, an exact scaling equation for the third order velocity correlations entering in 
the helicity flux definition has been derived under two hypothesis: (i) there exists a non-vanishing 
helicity flux, (ii) the flux becomes Reynolds independent in the limit of FDT. This relation predicts a 
r 2 -scaling for the particular third order velocity correlation entering in the definition of helicity flux, 
at difference from the celebrated linear behavior in r showed by the third order velocity correlations 
entering in the definition of energy flux. 

This simple fact tells us that different velocity correlation with the same physical dimension but 
with different tensorial structure may show different leading scaling properties. 

Moreover, even if overwhelming evidences indicate that the main physics is driven by the energy 
transfer, there can be some sub-leading new intermittent statistics hidden in the helicity flux prop- 
erties. 

Homogeneous and isotropic turbulence has, by definition, always a vanishing mean helical flux. 
Nevertheless, both fluctuations about the zero-mean, in isotropic cases, and/or net non-zero fluxes, 
in cases where inversion symmetry is explicitly broken, can be of some interest for the understanding 
of fully developed turbulence. 

In this letter, we analyze the helical transfer mechanism in dynamical models of turbulence 
HTlj H, [l^], built such as to explicitly consider helicity conservation in the inviscid limit. 

In a previous publication || we found the first strong numerical evidence that a Reynolds- 
independent helicity flux is present in cases where the forcing mechanism explicitly breaks inversion 
symmetry. In this paper we present a phenomenological argument which supports the idea that 
helicity behaves as a quantity passively transferred toward small scales by the mean energy flux, in 
agreement with both numerical findings in the true Navier-Stokes eqs. [|7| and with our previous 
numerical evidences about the strong intermittent properties of the helicity flux ||. 

In the following, we briefly summarize the main motivation behind the introduction of Shell Mod- 
els for turbulence. We summarize the main results obtained in || and we present the new argument 
about the passive character of helicity fluctuations. 



2 Shell Models 

Shell models have demonstrated to be very useful for the understanding of many properties connected 
to turbulent flows [jnj-[]nj|. The most popular shell model, the Gledzer-Ohkitani-Yamada (GOY) 



model ([0-[|T^]), has been shown to predict scaling properties for ((p) (for a suitable choice of the 
free parameters) similar to what is found experimentally. 

The GOY model can be seen as a severe truncation of the Navier-Stokes equations: it retains 
only one complex mode u n as a representative of all Fourier modes in the shell of wave numbers k 
between k n = k 2 n and k n+ i. 

It has been pointed out that GOY model conserves in the inviscid, unforced limit two quadratic 
quantities. The first quantity is the energy, while the second is the equivalent of helicity in 3D 



turbulence |Tj|. In two recent works || [12fl the GOY model has been generalized in terms of shell 
variable, u+,u~, transporting positive and negative helicity, respectively. It is easy to realize that 
only 4 independent classes of models can be derived such as to preserve the same helical structure 
of Navier-Stokes equations ||. All these models have at least one inviscid invariants non-positive 
defined which is similar to the 3D Navier-Stokes helicity. In the following, we will focus on the 



intermittent properties of one of them which has already been extensively investigated (see [12] , p20 
for more details). The time evolution for positive- helicity shells reads | jP2"| : 



w+ = ik n [An [u, u])* - vk 2 n u\ + 5 n>n J + , (3) 



with the equivalent eqs, but with all hehcity signs reversed, tor u . In (|d|J v is the molecular viscosity, 
/ + , /~ are the large scale forcing and A[u, u] refers to the non-linear terms of the model. Namely: 

A n [u, u\ = u~ +2 w++i + hu~ +l u+_ x + czu~_ x u~_ 2 . (4) 

— 5/12, C3 = —1/24 there exists two global inviscid 



It is easy to verify that for the choice 63 
invariants [|E|: the energy, 
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£ = E(k + l 2 + 



and helicity, 
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The equations for the fluxes throughout shell number n are 
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where Ei and Hi are the energy and helicity of the i-th shell, respectively: Ei 
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. Ei n and i/j n are the input of energy and helicity due to forcing effects, 

E m = »((/+(«?)• + /-(%)*)), H m = (/+«)* - /-(%)*)). 

In (^|) and ([]) we have introduced the triple correlation: 



uuu 



uuu \ 



(A+ +1 + A; +1 ) + (6 3 + l/2)(A ) t + A r :; 
(A+ +1 -A; +1 ) + (6 3 + l/4)(A+-A- 
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(7) 
(8) 

(9) 



Assuming that there exists a stationary state we have 



-T1 H 



0, where IT^ = k n ((uuu)^) 

and 11^ = k^({uuu)n) . Moreover, in the inertial range we can neglect the viscous contribution in 
(||) and (0), obtaining: 



In 



(uuu);) = k~ l E in , ((uuu);) = k~ 2 H in . (10) 
we have shown that in the case of a large scale forcing breaking inversion symmetry (i.e. 



/ + 7^ /~) energy and helicity fluxes coexist in the systems, both being Reynolds-independent. In 
the inertial range we may therefore neglect the viscous contribution and we obtain for the energy- 
triple-correlations and for the helicity-triple-correlation: 



[uuu)®) ~ k n \ 



(uuu);) ~ k n 2 . 



Ill 



Relation (|T0l) is the equivalent of what found for helical Navier-Stokes turbulence in |3|, [T(J . 

Let us remark that the coexistence of both energy and helicity fluxes is only possible due to the 
non-positiveness of helicity; in 2D turbulence, for example, a similar result, concerning enstrophy 
and energy cascades, is clearly apriori forbidden. 

In @ we have measured the statistics of energy and helicity transfers, by defining the two sets of 
scaling exponents: 



,(p) 



UUU) 



(p/3)\ 



it 



(12) 
(13) 



0.5 - 







4 



6 



8 



10 



Figure 1: Anomalous exponents for the helicity flux, 
1 9 obtained with N = 26 and v = 2 • 1(T 9 . 
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As one can see in Figure 1 we found that the even part of the two statistics coincides, i.e. 
C(2p) = ip(2p). On the other hand, the scaling exponents of odd moments are different. In particular 
the helicity exponents show a strong intermittent shape reminding to what one finds, for example, 
for the statistics of a passive scalar advected by a turbulent flow. Phenomenological and numerical 
evidences supporting the idea that helicity may be safely be considered as a passive advected quantity, 
at small enough scales, have already been presented in 0. 



3 Helicity advection 

Let us push further this idea by analyzing in more details the equation of motion describing the 
evolution of shell variables. In order to highlight the energy and helicity dynamics it is better to 
project the equation of motion into two new variables feeling the even and the odd part of the 
statistics respectively: 

w n = ^^ \ n = ^^. (14) 

Let us notice that due to the transformation properties with respect to the symmetry w+ <-> u~ it is 
obvious that w n will mainly feel the energy statistics while X n will be highly sensible to the helical 
properties of the flow. The total energy and helicity in the shell n will be E n oc \w n \ 2 + |A n | 2 and 
H n oc U(w n X* n ). 

Writing now the inertial-time evolution for these two variables we obtain: 

Wn = {w n +2W n +l + hWn+iWn^x + C 3 W n _iW n _ 2 ) + (15) 

( — A„+2A„+i — &3A n +iA n _i + C3A n _iA„_2j 



A n = («W2A n+ i + 63«WiA„_i - c 3 «v_iA n _2) + 

( — A n+2 1t; n+ i — &3A n+ iU>„_i — C3A n _iU>„_2) • 



(16) 



Equations ([TJ) and (|T6|) are identical to the set of coupled equations describing - in the shell world 
- an active scalar A n advected by a turbulent velocity field w n . The active character of A n is clearly 



dictated by the quadratic terms m the KHb 01 (|15|J . 



We now want to argue that at small enough scale the 0(A 2 ) part in the RHS of (|15|) becomes 
negligible and therefore we end up with a set of eqs describing the passive advection of the field \ n 
by the velocity field w n . 



Let us start by noticing that for the variables A 1 * 1 which enter in the definition of the fluxes we 
have: 

(17) 

and 



A + + A oc w n+ i [w n w n -i - A n A n _i] + A n+ i [iu n A n _i - A n w n _i] 
A + - A~ oc w n+ i [w n A n _i - A n w n _i] + A n+ i [w n w n _i - A n A n _i] . 



Furthermore, by imposing now two general scaling laws for the field w n ~ k n a and A 
obtain: 

A+ + A" ~ wl + w n Xl + w 2 n X n ~ k~ 3a + k~ a ~ 2p + k~ 2a - p 



„ ■ - kj we 



'uuuf n > ~ 



and 



A+ - A- ~ w 2 n X n + w n \ 2 n ~ k~ 2a ~P + k- a ~ 2 ?. 



(19) 
(20) 



By requiring now the asymptotic matching of fll9|) and of ( p0|) with the known behavior for the 
energy and helicity fluxes we have: 



( \i 



In order to satisfy both matchings we must to require that 

w n ~ A;" 1 / 3 A n ~ k~ A '\ 
1 



(21) 
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(23) 
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Figure 2: Log-log plot of Energy, E(k), and helicity, H(k), spectra compared with theoretical expec- 
tations. 



By plugging the above scalings in the equation of motion we discover that the time evolution for 
the field w n is mainly governed by the field itself, i.e. the coupling with the field A„ is sub-dominant 



at small scales and therefore as a direct consequence that the held \ n is passively advected by w n . Ui 
course the phenomenology here discussed will be affected by intermittency, i.e. the above exponents 
may show some weak/strong deviations from the predicted values, without, however having the 
possibility to change the role of dominant and sub-dominant contribution in expression ([H]) and 

OH)- 

In order to test this prediction we plot in Figure 2 the energy and helicity spectrum respectively 
E{k n )/k n ~ w\ ~ k~ 2 / 3 and H{k n )/k n ~ w n \ n ~ k~ 5 ^ 3 . As one can see, a part from sthe strong 
bottleneck effect showed by the helicity spectrum, the agreement is perfect. 

In conclusion we have investigated the Helicity transfer statistics in a shell model of turbu- 
lence. In particular, we have found that the strong intermittent properties of Helicity flux -measured 
numerically- can be explained in terms of the phenomenological evidence that Helicity becomes a 
passive quantity advected -at small enough scale- by the energy flux. 
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